ACT-03/01, CTP-TAMU-11/01, |Ei^h/010334S 
Inflationary Baryogenesis 



RAGHAVAN RANGARAJAN^") AND D.V. NANOPOULOS^'')'^'^)'^'^) 
Theoretical Physics Division, Physical Research Laboratory 
Navrangpura, Ahmedabad 380 009, India 
Department of Physics, Texas A&M University 
College Station, TX 77843-4242, USA 
Astroparticle Physics Group, Houston Advanced Research Center (HARC) 
Mitchell Campus, The Woodlands, TX 77381, USA 
^'^^ Chair of Theoretical Physics, Academy of Athens 
Division of Natural Sciences, 28 Panepistimiou Avenue 
Athens 10679, Greece 

Abstract 

In this letter we explore the possibility of creating the baryon asymmetry 
of the universe during inflation and reheating due to the decay of a field asso- 
ciated with the inflaton. CP violation is attained by assuming that this field is 
complex with a phase that varies as the inflaton evolves. We consider chaotic 
and natural inflation scenarios. In the former case, the complex decaying field 
is the inflaton itself and, in the latter case, the phase of the complex field is 
the inflaton. We calculate the asymmetry produced using the Bogolyubov 
formalism that relates annihilation and creation operators at late time to the 



annihilation and creation operators at early time. 



Typeset using REVT^ 



I. INTRODUCTION 



Explaining the origin of the matter-antimatter asymmetry of the universe is an essential 
ingredient in our understanding of the history of the universe. In this article, we study the 
possibility of creating the baryon asymmetry of the universe by the production of particles 
during inflation and reheating by the decay of a complex field related to the infiaton. We 
consider the case where the complex decaying field is the infiaton itself as well as the case 
where the phase of the complex field is the infiaton, as in natural inflation. The former case 
is similar to chaotic inflation but with a complex infiaton. By assigning baryon number to 
a scalar field present during inflation and introducing a baryon number violating coupling 
between this field and the infiaton we find that there is a net baryon number asymmetry in 
the produced particles. 

The notion that the infiaton plays a role in baryogenesis is not new. If the reheat 
temperature is above the mass of certain heavy particles, such as GUT gauge and Higgs 
bosons, then the latter are thermally produced and their subsequent out-of-equilibrium 
decays create a baryon asymmetry. The production of heavy GUT gauge and Higgs 
bosons or squarks by the direct decay of the infiaton when the reheat temperature and/or 
the infiaton mass is less than that of the heavy bosons has also been considered. Once again, 
the out-of-equilibrium decays or annihilations of these particles gives rise to the baryon 
asymmetry of the universe. In all the above scenarios CP violation enters into the 

couplings of the heavier bosons to lighter particles. In our scenario the baryon asymmetry 
is produced in the direct decay of a field associated with the infiaton. Furthermore the CP 
violation must manifest itself in the decay of this field. In this respect it is similar to the 
scenario mentioned in Ref. and discussed in more detail in Ref. in which the baryon 
asymmetry is created by the direct decay of the infiaton. However, unlike in Ref. 0, in 
our scenario CP violation is provided dynamically through the time dependent phase of an 
evolving complex infiaton or of a complex field associated with the infiaton. We explicitly 
calculate the asjTiimetry in our scenario and compare it to the baryon asymmetry of the 
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universe. We follow the work of Ref. [§] in which the asymmetry was calculated in the context 
of a universe that contracts to a minimum size, bounces back and then expands. The universe 
was static at both initial and late times. The B- violating coupling was XR{(f)*Ailj + ?/'*A*0), 
where R is the Ricci scalar, (p carried baryon number and A was a complex function of time, 
which provided the necessary CP violation. In this work, we have adapted the formalism of 
Ref. to consider the asymmetry that might be created in a more realistic universe that 
is initially inflating and then enters a reheating phase followed by the standard evolution of 
the universe. Furthermore we have given a more realistic source of CP violation, namely, a 
time varying complex field. 

Recently Ref. appeared in which the authors discuss the generation of the baryon 
asymmetry during preheating in a scenario similar to the one discussed here. We discuss 
later the differences and similarities between our work and theirs. 

As in Ref. 0, we consider a lagrangian consisting of two complex scalar fields and 
ip. (f) and %l) are assumed to carry baryon number +1 and respectively and we assume a 
B- violating term 



where A is a dimensionless constant and rj is related to the infiaton field and is complex. The 
baryon number of (f) and if) particles is established by their interactions with other particles 
in the Standard Model. Q The latter are not included in our lagrangian below as they do 
not enter into our calculations. We assume that the initial velocity of the rj field and/or the 
shape of its potential ensures that its phase varies as the infiaton rolls down its potential. 
Thus we have dynamic CP violation. 

To obtain the asymmetry in our scenario we use the fact that the annihilation and 



To ensure that the baryon asymmetry created is not erased by sphaleron processes, we assume 
that the interaction in Eq. (|l]) also violates B-L. This may be achieved, for example, if ip carries 
no lepton number. 




(1) 
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creation operators for the fields and ijj are not tlie same during tlie inflationary phase and 
at late times after reheating. However, the annihilation and creation operators at late times 
can be written as linear combinations of the annihilation and creation operators during the 
inflationary phase using the Bogolyubov coefficients. So, 

4 = A,at + Bkbti + A' at + B',h% , (2) 
= C^a\ + ^fc&f + C'a\ + D'.ht^ , (3) 

where and 6^ are operators at late times and and 6^ are operators at early times in the 
inflationary phase. Similar expressions exist for and 6^. In the Heisenberg picture, if we 
choose the state to be the initial vacuum state, then it will remain in that state during its 
subsequent evolution. One can then see that the number of particles and antiparticles of 
momentum k, given by (0|a^^a^|0) and (0|&k^6k|0) respectively, are non-zero and proportional 
to + |5^|2 and |Cfc|2 + Furthermore, if \Bk\'^ + ^ jCfep + jC^p then one 

obtains a baryon number asymmetry. (If A = in Eq. |I]then B' and C are and \B\ = \C\, 
and one gets no asymmetry.) 

Particle number, and correspondingly annihilation and creation operators, are well de- 
fined only in adiabatic vacuum states. However, the vacuum state in the inflationary era, 
i.e., our in state, must be chosen judiciously to avoid infrared divergences. This is discussed 
in Section HI. 

The framework of this article is as follows. In Section II we present the lagrangian density 
for the complex scalar fields and ip relevant to our calculation and obtain their equations 
of motion. We then write down the Fourier decomposition of and ip during the inflationary 
phase and during reheating. General expressions for the coefficients A — D' in Eqs. ^ and 
^ have been derived in Ref . |^ . We shall present these results without rederiving them and 
then present the general result for the baryon asymmetry of the universe. In Section HI, 
we present the particular solutions for our scenario of a universe that undergoes exponential 
inflation (a ~ e ) followed by an inflaton-oscillation dominated phase (a ~ 

t^/S). We then 

calculate the total baryon asymmetry for this scenario in the context of chaotic and natural 
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inflation. We conclude in the last section. In the Appendix we discuss issues related to the 
regularisation of infrared divergences and the necessity of an infrared cutoff to satisfy the 
conditions of perturbation theory. 

II. 

Consider a lagrangian density 

L = ^/=^[^?^'^5^0*5,0 + g'^'d^rO,^ - {ml + i^RW<P - (m^ + i^R)^^^ (4) 
+g^''d^r)*d^r] - m^r]*r] - V{r]) - X{rf(t)*'il: + r?*V»], (5) 

where m<^_y, are the masses of the respective fields and ^<^_^ are their couplings to the cur- 
vature. We have assumed that r] is minimally coupled. ¥{7]) includes all interactions of r] 
other than the coupling to (p and ip already listed above. Below we shall consider natural 
infiation and chaotic infiation scenarios. In the former case the inflaton will be associated 
with the phase of r^. In the latter case we assume that the complex r] field is the infiaton 
field. (Thus we are really considering an extension of chaotic inflation since the inflaton is 
now complex.) We shall assume a spatially flat Robertson- Walker metric. The equations of 
motion for the above flelds are 

+ 3(-)0 - 4v'0 + K + ^^^)'^ + VV' = 0, (6) 
V; + 3(-)^ - ^V^V^ + {ml + i^R)'iP + \ri*^ct> = 0. (7) 

We now write 

with k = ^{nx,ny,nz) and a similar expression for ip{x). The equations of motion for the 
Fourier coefficients 0k and ■^k are (note that 0k and ■^k are operators) 

Mt) + (K-j-2-l- + ^l + ^'^R)Mt) + XR^Mt) = 0, (9) 
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Mt) + + + e^^)^k(t) + XRA*Mt) = 0, (10) 



0k and ■j/'k satisfy 

[0k(t),0k'W] =^5k,k', (11) 

[Mt),ri.'{t)] = ts^y. (12) 

To use the results derived in Ref. ||^ we shall have to solve Eqs. |^ and |T0| for times 
during inflation and during reheating. To simplify our calculations we shall assume that the 
fields (p and ip are massless and minimally coupled, i.e., m^^^ = 0,^^,^ = 0. (Typically, a 
spin zero particle will obtain a mass of order H during inflation, if H is greater than its bare 
10| . However, this does not occur if the mass is protected by a symmetry.) 



mass 



To facilitate the use of perturbation theory and to be able to define an in state we assume 
that the B- violating interaction switches on at some time ti. Let t2 be the time when inflation 
ends, ts be the time when reheating ends and t/ be the final time at which we evaluate the 
baryon asymmetry. We assume that B- violation vanishes after ^3. The annihilation and 
creation operators aX tf can be expressed as linear combinations of the annihilation and 
creation operators at an early time ti before ti in the inflationary era. These relations have 
been derived perturbatively to order in Ref. giving 

a?,k = [4a + ^>^'Ht) - iX'PtHt]ai^ + (1 - iX^Hf) + iX'4Ht]l^U 

- zX[ath - f3th]at^ - ^X[ath - f^thW,-^ , (13) 
bt = [«f(l - ^^'Ht) + ^X'^3tHt]btl + [/?,^(1 + zX^Hf) - zX'atH^]al^ 

+ ^X[ath - Pthjbfl + ^X[ath - PthjaU , (14) 
4,k = [«t(l + ^>^'Ht) - iX'PtHt]at^ + [/?r (1 - iX'Ht) + iX'atHt]ht^ 

- - ^tn]aU - ^>^[4i3 - PtmbfU ' (is) 

bt^ = [«r (1 - ^X'Ht) + iX^PtHtm + [/3,^(1 + iX^Ht) - ^X'atHt]aU 

+ ^A[ar/: - Pti;]btl + ^X[atI; - f3th*]aU , (16) 



where 



h= iZodtv'mtmtit) , (17) 

h= IZ.dtv\t)xt{t)xt{t), (18) 

h= I-o.dtv\t)xt{t)xt{t) , (19) 

h= IZodtv'{t)xt{t)xt{t) , (20) 



and 



Ht= IZodtJZodt'v'mtit) 

X Atit,t')xtit')v*\t') , (21) 

H^= IZodtJZodt'v'mtit) 

x/\tit,t')xt{t'W\t') , (22) 

Ht= IZodtJZodt'v'mtit) 

X ^tit,t')xt{t'W'{t') , (23) 

Ht= IZ.dtJZ.dt'^\t)xt{t) 

x^t{t,t')xt{t')ri*\t') , (24) 

(25) 

and the Hf are defined as Hf with replaced by r/*^ and , x^* and A"^ replaced by xt, 
xt* and A*^ respectively. Above, xt and xt are complex functions (not operators) that solve 
Eqs. ^ and [l^, with A set to 0, respectively. A^ and A^ are the retarded Green's functions 
for Eqs. ^ and |10| respectively, i.e., they satisfy Eqs. ^ and ^ with A set to and a delta 
function 6{t — t') on the r.h.s. of the equations. The subscripts on the coefficients and 
(3k and on the functions Xk and A^ refer to |k|. and (3^ are complex and satisfy 

\4?-\Pt? = \4?-Wt? = ^- (26) 

We assume that the initial state at ti is the vacuum state. Then, in the Heisenberg 
picture, the number of (j) particles of momentum k at is given by 



= \(it? + >^'[\4?\h? + \Pt\\\h?) + 2Re{atpt{^Ht - hit))] (27) 

and 

= \Pt? + I VsP + l/?.1'(|/2n + 2Re{at(3'^-\^Ht + AVa))] (28) 

where we have used 

2/mi7i = - |/2|^ (29) 

2ImHi = \h\^ -\h\^ . (30) 

The baryon asymmetry for particles of momentum k at is thus 

/\Nt{tf) = {N^itf)) - {N^itf)) = X\\h\' - \h\') , (31) 



where we have used Eq. ^ and 

- Ht = iIJl - ihll . (32) 

The reader is referred to Ref. p[ for a more detailed derivation of the above results. Note 
that the asymmetry does not depend on at and (3k implying that the asymmetry is indepen- 
dent of the purely gravitational production of particles in the expanding universe. (^ = 
does not imply a conformally invariant universe. Therefore there is non-zero purely gravita- 
tional production of particles in our scenario but it does not contribute to the asymmetry.) 

To obtain the net baryon number at we sum over all momentum modes and take the 
continuum limit. Since we ultimately wish to express the baryon asymmetry as the baryon 
number density to entropy density ratio, and the baryon number does not change after ^3, 
we write 



^Eqs. (4.1c) and (4.9) of Ref. Q contain typographical errors. The corrected equations are dis- 
played here in Eqs. |l5| and 32 respectively. 
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nsih) = lim^^oo l/{[Laits)Y) E^i^N^) 

= l/{2n^[a{h)f) /o°° dkk\/\Nt) . (33) 

Working in the approximation that at ^3 the inflaton completely decays and the universe 
instantaneously reheats to a temperature T3, the baryon asymmetry of the universe is 

\2 /■oo 

BAU ^ ns/s = / dkk\\h\^ - |/3n/[(27rV45)^7*T|) , (34) 

where we have assumed that there is no dilution of the baryon asymmetry due to entropy 
production during the subsequent evolution of the universe. Note that because the effective 
coupling is a complex function of time the baryon asymmetry is obtained at O(A^) and not 

For standard reheating, ^3 F"^, where F is the dominant perturbative decay rate of 

2 

the inflaton. We take F = |^mj„j, corresponding to the decay of the inflaton of mass minf 

I \- 1 — - 1 

to some light fermion-antifermion pair. Furthermore, T3 = y-^^y pi^s)^ = O.Gg^, *(Mp;F)2 

TT| , where p here refers to the inflaton energy density. pits) ~ p{t2)[a(t2) / a(ts)]^ ■ We 



assume that the reheat temperature is not high enough for GUT B- violating interactions to 
be in equilibrium and wipe out the asymmetry generated in our scenario. On the other hand, 
we do not restrict ourselves to reheat temperatures below 10^ GeV to avoid the gravitino 



problem [113] as we consider the possibility that the gravitino might be very light. 



III. 

To obtain the baryon asymmetry, we need to evaluate I2 and I3. This requires obtaining 
xt and xt- shall need to perform the integral for I2 and I3 only from ti to ^3 as B- 



The final temperature T4 at the end of reheating is also a function of the interactions of the 



inflaton decay products which we have ignored. See Ref. |12| and references therein for a discussion 
of thermalisation of the decay products. 
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violation vanishes earlier than ti and after ^3. Solutions of Eqs. ^ and ^ for A = and 
a(t) = a-f^, (c 7^ 1,-1/3) and a{t) = ae^* have been obtained in Ref. [0. Using them we 

get g 



Xt"" = [a{t)] 



Cl( 



3a{ty^H 



a{t)H' 



+ C2( 



-1 



?,a{tfH 



Xk 



,1/2^(2) 



a{t)H' 
-3kt. 



^^\{tf' f ^ a{t) ' ^^a{tf' "§ ' a(t) 
Hy{z) are Hankel functions. The commutation relations Eqs. |Tl| and ^ imply that 



for ti<t <t2, (35) 
for t2<t<h. (36) 



Xk Xk - Xk Xk = ^ 



(37) 



and therefore the constants Ci^2 and c[ 2 satisfy 



I |2 I |2 
C2 - Ci 



I ' |2 _ I ' |2 

1^2! |Ci| 



= 37r/4, 
-37r/4. 



(38) 
(39) 



The constants Ci and C2 define an initial vacuum state in the inflationary era. If one 
makes a choice of the de Sitter invariant vacuum state (ci = and C2 = ^JSn/A) as the 
in state then it is well known that such a state suffers from an infrared divergence. One 
option then is to choose the constants ci and C2 appropriately so as to cancel the infrared 



divergences even though such states will no longer be de Sitter invariant. In Ref. |14| it is 
suggested that one may choose the constants Ci and C2 as below so as to cancel the infrared 
divergences: Q 



C2 



c, = {kt2/a{t2))-r 
{{kt,/a{t,)r^^A'l\ 



(40) 
(41) 



^ Note that the definition of the mode functions Xk in Ref. |8| differs from that in Ref. [|T4| by a 
factor of a^/^. 

^ Other mechanisms have also been suggested to avoid the infrared divergences. For example, see 
Ref. 115 1 . 
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with p > 0. (We point out in the Appendix that there is also an upper hmit on p that 



was not mentioned in Ref. [|14[-) As we discuss in the Appendix, the nature of the infrared 
divergences is shghtly different in our case. Though the above choice for the constants ci^2 
with appropriately chosen values of p would make the final integral over k infrared finite, 
the integrands for the intermediate integrals over t become very large for small values of k, 
irrespective of the value of p. This leads to a problem with perturbation theory since the 
latter requires that A^|/2,3p should be less than 1. ^ This is discussed in more detail in the 
Appendix. Therefore we are forced to introduce a low momentum cutoff ki to justify our 
use of perturbation theory. We choose ki such that kL/a2 = 1/^2- Since the low momentum 
cutoff automatically regulates the integral over k, we then choose Ci = and C2 = ^JSn/A. 
Continuity conditions for 0(x) and a{t) at ^2 imply that x(t) and (i/(it(x(t)/[a(t)]^/^) 

are continuous at ^2, and these boundary conditions then give us c[ and c'2. We have verified 
that the values of c[ and obtained from the continuity conditions satisfy Eq. ^ 

At this stage we need to specify i]{t). If we write r](t) as -^cr(t)e*^*^*^ (where a{t) is real), 
it is the time varying phase of 77 that provides the CP violation necessary for creating a net 
baryon asymmetry. 

Chaotic Inflation 

We first consider the case of chaotic inflation in which the rj field represents a complex 
infiaton field. In the absence of any potential for 6, the equation of motion for 6 is 

a^e + SHa'^e + 2&ae = . (42) 

In a more realistic model V{ri) will imply a potential for 6. The equation of motion for cr is 

a + 3Ha + rn^a - ^ V = . (43) 

We assume that Qif) evolves starting from an initial value of at t = tj and an initial 
velocity Qi. We choose Qi consistent with a universe dominated by the potential energy of o. 



^ We thank D. Lyth for pointing this out to us. 
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Therefore we take 9i = m/2. During inflation, <^ m^a^ and a ~ Mpi. Hence a/a ^ H 
and we ignore the last term in the equation of motion for 6. Then 

^^(t) = ^(l-e-3^(*-*»)) U<t<h (44) 

We take H to be constant during inflation and corresponding to the initial energy density 
of the universe with a(ti) = 3Mp;. From above, one can see that 6 is much less than m for 
most of the inflationary era and so we ignore the last term of Eq. ^ for this era. Invoking 
the slow roll approximation one may also ignore a during inflation. 

During reheating the a and the 6 fields are coupled and we can not ignore the last terms 
of their respective equations of motion. However, to obtain 6{t) it is simpler to first rewrite 
T] as ■^{'^1 + 11^2) ■ Then the problem reduces to one of two uncoupled damped harmonic 
oscillators with solutions, ni = {Ai/t)cos{mt + a) and k,2 = {A2/t)cos{mt + j3). Here we 
have assumed H = 2/{3t) during reheating. 6{t) is then tan'^{K,2/ i^i)- The constants Ai, 
A2, a and (3 are determined by the values of K12 and their time derivatives at ^2 which can 
be obtained from the values of 6 and a and their time derivatives at t2. We take t2 ~ 2/m, 
as the inflaton starts oscillating when 3H ^ m. a{t2) ~ Mpi/6. 

Eq. ^ implies that 6 becomes nearly constant within a few e- foldings after tj. If the B- 
violating interaction switches on subsequent to this then 6 is approximately constant between 
ti and ^2- Furthermore, since ^(^2) is practically zero, there is practically no rotational motion 
during reheating in the absence of any potential for 6. So during reheating 6 takes values 
of 62 and $2 + TT during different phases of the oscillation of a, where 62 is the value at t2- 
{6 changes discontinuously at the bottom of the potential where a is 0.) Since the relevant 
phase in I2 and I3 is 26 the above implies that the CP phase is practically the same for the 
interval ti to and hence one should expect very little asymmetry. 

Natural Inflation 

We now consider natural inflation, in which case a{t) = f where / is the scale of spon- 
taneous symmetry breaking in the natural inflation scenario. In the presence of an explicit 
symmetry breaking term that gives mass nig to the inflaton 6 the equation of motion for 6 
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IS 



fe + sHfe + mlfe = o . 



(45) 



We assume that is constant during inflation between ti and ^2 and is of 0(1). Our results 
are insensitive to ti for ti earlier than about 10 e-foldings before the end of inflation. At 
t2 ~ 2 /me when 3if ^ m the 9 fields starts oscillating in its potential. Between t2 and 9 
evolves as 



Now we obtain the asymmetry numerically. Smaller the value of longer is the period of 
reheating contributing to the asymmetry. But for g < 10~^, I2 and become independent 
of g and then the g dependence in BAU enters through a{tz) and the reheat temperature 

T3. For g < IQ-^, we get 



As we have mentioned before, perturbation theory requires that A|/2,3p must be less than 
1. This translates into an upper bound on A of 10~^^. Then even for g = 10~^ we get 
insufficient asymmetry. Other values of g give even less asymmetry. 



In conclusion, we have discussed a mechanism for creating a baryon asymmetry dur- 
ing inflation and reheating. While the scenario illustrated above does not create sufficient 
asymmetry, it may be easily modified to accommodate a potential for 9 which can give rise 
to a much larger asymmetry. A possible potential for 6 for the chaotic inflation scenario 
is W{9) — m^(7^(l — cos^), which is equivalent to tilting the inflaton potential. Unlike in 
the analogous axion and natural inflation models, here both a and 9 would be varying with 
time. Hence such a potential may allow for chaotic orbits and so would have to be studied 
with care. 



9{t) = 9{t2)^cos[mg{t-t2)]. 

V 



(46) 



BAU = X^g{2 X 10^°) . 



(47) 



IV. CONCLUSION 
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We point out here that we include both the inflationary phase and the reheating phase 
in our calculation. Contributions during both phases do get mixed up in the evaluation 
of the asymmetry because of the presence of and I/3P, where the time integrals in I2 
and J3 include both the inflationary and the reheating eras. In fact we find in the natural 
inflation case that though the phase is taken to be constant during the inflationary era, the 
net baryon asymmetry for a fixed value of A decreases if we do not include the inflationary 
era in the integrals I2 and 1^. W This indicates that one should not ignore the inflationary 
era when calculating the asymmetry. 

While we were writing up this paper Ref. appeared. In this paper the authors discuss 
the generation of baryon asymmetry during reheating in a scenario similar to ours. The 
two calculations have some differences however. Our calculation is carried out in curved 
spacetime while the mode functions in Ref. are obtained in Minkowski space. We consider 
standard reheating while they consider the more complicated preheating scenario. In both 
calculations the source of CP violation is a time varying phase. The authors of Ref. 
suggest that the CP violating potential for the baryonic fields may be induced by their 
direct coupling to the inflaton or through loop effects involving the baryonic fields and other 
fields and then presume a form for the phase. We provide a specific scenario in which the 
infiaton, or a field related to the infiaton, which is coupled to the baryonic fields, is complex 
and its time varying phase dynamically provides CP violation. Involving the infiaton and 
its phase seems to us to be a simple and a very natural approach to obtain a time dependent 
phase. Our calculation includes both the infiationary and the reheating eras which, as we 
have pointed out, seems to be appropriate for our case. 



^The values of |/2,3p also decrease and so, in principle, one may obtain a slightly larger asymmetry 
by choosing a larger value of A which is still consistent with the use of perturbation theory. 
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APPENDIX 



In Ref. [14| infrared divergence in for Robertson- Walker universes with power 

law expansion and exponential expansion is discussed. (In this Appendix, is a generic scalar 
field.) Now 



where Xk^i^) mode functions as defined in Ref. (they are a^/^ times the mode 
functions defined in this paper) and t = /* a^^{t")dt" . If one argues that H^l\z) ~ Hf\z) ~ 

2 2 

iYi(z) ~ z~^^'^ as 2; ^ the integrand in the above integral ~ |ci — C2pA;"'^, which implies 
a divergent 2-point function unless one chooses Ci{k) and C2{k) such that |ci — C2I -^0 fast 
enough as k 0. Hence the choice given in Eqs. (|40| ) and (^) is suggested in Ref. |jT4 



which eliminates the infrared divergences if p > 0. We would like to point out here that 
there is also an upper bound on p. Unless p is less than 3, the term proportional to Jf in 
the 2-point function, namely, ~ |ci + C2pJi ~ will also diverge at low values of k. 

Unlike in Ref. [1^ the divergent quantities for us are ~ / d^kxt*xt* ~ / d'^f^xtxt 



and the form of the divergences is slightly different. Adapting the arguments of Ref. |T4 



we would choose | < p < | for the constants Ci and C2 given in Eqs. ( |iO| ) and (^TJ) to avoid 
infrared divergences in the integral over momentum 

In our case we also have intermediate integrals I2 and I3 over time and this leads to an 
additional problem. The k dependence in the integrand of, for example, I3, is contained in 
— (ci — C2YY^{z) + (ci + C2)^Jf{z) + 2i{c\ — C2)Ji{z)Yi{z), where z is the relevant argument 
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of the Bessel functions depending on whether the universe is in the inflationary or the 
reheating era. The k dependence for I2 is similar. When k is small the argument of the 
Bessel functions becomes small. For low k values the first term goes as k^P~^, the second 
as k^~'^P and the third is k independent. Now perturbation theory requires that A^|/2,3p be 
less than 1. However, since 2p — 3 and 3 — 2p cannot both be greater than 0, I2 and Is will 
become very large at low k values. Thus, without a low momentum cutoff, perturbation 
theory breaks down at some point for any finite value of A. We emphasise again that this is 
an issue related to the validity of perturbation theory and not to the infrared divergence of 
the integral over momentum. The latter can be regulated by the choice of constants ci and 
C2 mentioned above, irrespective of whether or not perturbation theory is valid. 
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